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Abstract 

Since the 1970's, theories of deformation and failure of amorphous, sohdUke materials have 
started with models in which stress-driven, molecular rearrangements occur at localized flow defects 
via "shear transformations." This picture is the basis for the modern theory of "shear transforma- 
tion zones" (STZ's), which is the focus of this review. We begin by describing the structure of the 
theory in general terms and by showing several applications, specifically: interpretation of stress- 
strain measurements for a bulk metallic glass, analysis of numerical simulations of shear banding, 
and the use of the STZ equations of motion in free-boundary calculations. In the second half of 
this article, we focus for simplicity on what we call an "athermal" model of amorphous plasticity, 
and use that model to illustrate how the STZ theory emerges within a systematic formulation of 
nonequilibrium thermodynamics. 
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I. OVERVIEW 



Deformation and failure of amorphous, solidlike materials is a large topic with a long 
history. This class of materials includes structural and metallic glasses, glassy polymers, 
dense colloidal suspensions, many kinds of granular materials, and a huge range of biological 
substances. Despite the challenge presented by their inherent disorder, these materials have 
received special attention because of the broad range of applications in which it is essential 
to predict their failure modes. The ways in which amorphous solids behave when driven 
by external forces have been studied throughout the latter half of the 20'th century by 
disparate groups of physicists and materials scientists. The field has advanced significantly 
in recent years, largely because of insights gained from numerical simulations and partly 
because of the introduction of a few new theoretical concepts; but even the basic form of an 
acceptable, first-principles theory has remained a matter of debate. 

Among the computation-inspired advances is the "shear-transformation-zone" (STZ) the- 
ory of amorphous plasticity, which we proposed in 1998 [2I, and which has gone through a 
series of modifications and extensions since then. (For example, see 0-9|.) The STZ theory 
is the principal topic of this review. We emphasize two main points: first, that the STZ 
theory starts with specific assumptions about the nature of molecular rearrangements to 
arrive at testable predictions of experimental data; and, second, that the theoretical hy- 
potheses are strengthened by being developed in the context of fundamental principles of 
nonequilibrium thermodynamics. 

From its inception ^ , the STZ theory was intended to be an extension of the flow-defect 



theories of TurnbuU, Cohen, Spaepen, Argon and others |10l-ll4j| in which localized clusters 
of molecules undergo irreversible rearrangements in response to applied shear stresses. In 
a very rough sense, the flow defects in amorphous materials play the role of dislocations in 
crystals by being the agents of plastic deformation. All of these theories, including STZ, 
start by assuming that the material of interest is solidlike - that it has a shear modulus - 
and that the flow defects or dislocations allow it to behave in some respects like a liquid. 



Thus, the fiow-de 



'ect theories are qualitatively different from fluid-based theories such as 



mode-coupling |l5l-ll7l|. which starts from a liquidlike, many-body description and predicts 



the onset of solidlike behavior at high densities and low temperatures. 

The most important way in which the STZ theory departs from its predecessors is by 
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recognizing that the flow defects must possess internal degrees of freedom. The STZ's not 
only appear and disappear during conflgurational fluctuations; they also transform from one 
orientation to another, and the net rate at which these transformations occur determines 
the rate of irreversible shear deformation. As we show later in this review, the equations 
of motion for this orientational degree of freedom predict that the system undergoes an 
exchange of dynamic stability from jammed to flowing states at a yield stress or, more 
accurately, at a minimum flow stress. This dynamic interpretation of yielding is not exactly 
the same as the conventional picture of stress-driven barrier crossing; it has as much to do 
with entropy generation as it does with deterministic dynamics. The orientational degree of 
freedom also gives the system a directional memory; for example, it produces Bauschinger 



effects. 



18| 



In recent years, the STZ theory has acquired a thermodynamic complexion with the 
introduction of an effective disorder temperature. Like the flow-defect picture, the effec- 



tive temperature idea has historical roots in Cohen and TurnbuU's theories of glasses 19 |. 
Other investigators, notably Spaepen , described the intrinsically disordered state of non- 
crystalline materials by a free volume vj. Those authors perceptively recognized that the 
relevant deflnition of is not as an extensive excess volume measured from some densely 
packed state, but as an intensive quantity - the inverse of the derivative of a conflguratio 
nal entropy (i.e. a dimensionless entropy associated with the mechanically stable molecular 
positions, without kinetic or vibrational contributions) with respect to the volume. Thus 
they proposed that the density of flow defects might be proportional to a Boltzmann-like 
factor, exp (— constant/f j), and not just to vj itself. This thermodynamic analogy was fur- 
;her developed for strictly athermal materials such as powders by Edwards and coworkers 
20, l2l|] . Lemaitre 22|, |23| has used the free volume in a way that is even more closely closely 
related to the effective-temperature analysis described here. 

Whenever the mechanical behavior of the system is determined by interactions between 
its elementary ingredients, it becomes important to couch the thermodynamic analysis in 
terms of energy rather than volume. Then the analog of the free volume f/ is an effective 
disorder temperature Te/-/, deflned here to be the derivative of the conflgurational energy 
with respect to the conflgurational entropy. T^ff equilibrates with the ambient temperature 
T at high T, but may fall out of equilibrium at low T when molecular disorder is generated 
by mechanical deformation, or when glass-forming liquids are quenched through a glass 
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transition. In the latter case, Tg// is the same as the Active temperature [24, |25| . In analogy 
to the free- volume formula, the steady-state STZ density is proportional to exp (— T^/Te//), 
where ez = ksTz is a characteristic STZ formation energy. 

Throughout this review, we consider only materials in which localized molecular rear- 
rangements occur in the presence of an ambient, high-frequency noise generated by either 
thermal or mechanical fluctuations. Thus, we do not consider strictly athermal situations 
that may occur in granular or colloidal materials very near the jamming transition. In par- 
ticular, the physics we discuss is quite different from that observed in athermal quasi-static 
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(AQS) numerical simulations 



3H33|. in which the system is relaxed to its nearest energy 
minimum after each of a sequence of small strain steps. Both AQS simulations, and fully 
dynamic simulations in the AQS limit of vanishingly small temperature and strain rate, 
characteristically exhibit system-spanning, crack-like events, and size- dependent spectra of 
stress fluctuations. 

Despite the obvious discrepancy between AQS and normal behavior, Lerner and Procaccia 



34j | have asserted that the system-spanning events observed at or near the AQS limit cast 
serious doubt on all plasticity theories that are based on the assumption of localized events. 
We always have held that such behavior may be typical of slowly driven, nearly jammed, 
granular materials; but conventional, amorphous molecular materials exhibit no such size 
dependent behavior. Since the models being simulated have normal, finite-ranged, molecular 
interactions, it is clear that there must be a crossover from AQS to normal behavior as a 
function of increasing temperature and strain rate. The system-spanning events must be 
suppressed by thermal and mechanical noise; and constitutive relations and fluctuation 



spectra must behave proper 
Procaccia and colleagues 



y in the thermodynamic limit of infinite volume. More recently, 
35| have demonstrated that this is indeed the case; in fact, they find 
that the expected crossover to localized events and normal plasticity occurs at exceedingly 
low temperatures. This observation makes it clear that it is inappropriate to overgeneralize 
the AQS results as in js^, and reaffirms that the existence of thermal or mechanical noise 
places amorphous molecular systems - to a very good first approximation - squarely in the 
regime where local events control dynamics, consistent with the STZ theory. 

A description of amorphous plasticity that we believe must be closely related to the 



STZ theory is known as "Soft Glassy Rheology" or "SGR" 



36 



37l |. This theory describes 



a broad range of glassy behaviors in terms of distributions of localized, noise activated. 
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displacement processes. Like the STZ theory, SGR exhibits transitions between different 
kinds of jammed and flowing states and can account, at least qualitatively, for a variety of 
phenomena observed in soft materials. However, SGR begins by assuming a fixed, broad 
distribution of activation energies and the existence of a noise temperature that controls the 
activation rates. The meaning and dynamics of the SGR noise temperature remain unclear 
at present; we do not even know whether it might be same as the STZ effective temperature. 
In contrast, the STZ theory is based on a specific model of molecular rearrangements that 
have been observed directly in numerical simulations and analog experiments. Both the 
effective temperature and the dissipative processes that generate it are defined in terms of 
those molecular degrees of freedom. 

This review consists of two main parts, presented in SecsJITjand lllli In Secjlll we exhibit 
one limiting form of the STZ equations of motion and very briefly describe several appli- 
cations, specifically: interpretation of stress-strain measurements for a bulk metallic glass, 
analysis of numerical simulations of shear banding, and the use of the STZ equations in free- 
boundary calculations. In Sec lIIIl we focus for simplicity on what we call an "athermal" 
model of amorphous plasticity, and use that model to illustrate the thermodynamic basis of 
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; but some 



the STZ theory. The thermodynamic arguments are taken primarily from 
features of this analysis have not appeared elsewhere. The derivations in Sec JIIII provide a 
first-principles rationale for the formulas used in SecJTTl 



II. SUMMARY AND SELECTED APPLICATIONS OF THE STZ THEORY 



A. Elasto-Plasticity 

Our first step in summarizing the STZ theory is to place it within a general set of 
Eulerian equations of motion for elasto-plastic deformation in a solidlike material. Consider 
a d dimensional system; let i, j... be spatial indices; and use summation convention. Then 
write the stress tensor aij in the form: 

= -p6ij + Sij, P= ^kk, (2.1) 

where p is the pressure and Sij is the traceless, symmetric, deviatoric stress. In analogy to 
fluid dynamics, let Vi{x,t) denote the material velocity at the physical position x = {xi} 
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and time t. The acceleration and continuity equations are: 

^ dvi_ _ doi^ _ _dp_ _^ dsij__ dpo d_ 

dt dxj dxi dxj ' dt dxi * 
Here, po is the density, and d/dt denotes the material time derivative acting on a scalar or 

a displacement or velocity field: 

d _ d d , , 

dt dt dxk 

For most purposes, we assume that po remains approximately constant and thus do not need 
the second, continuity equation in Eq. fl2.2p . 

The problem of disentangling elastic and plastic deformations has long been one of the 
more serious challenges in solid mechanics. This problem is not our main topic. (See [ssl for 
a recent review.) Rather than address it, we assume - as a simple working approximation 
- that the total rate-of-deformation tensor can be written as a linear superposition of 
elastic and plastic contributions: 

- 2 \dx^ + dx,) Vt \ 2K + 2/i )^ ^ ' 

where /i is the shear modulus, K is the inverse compressibility, and Dfj is the plastic rate 
of deformation. The symbol V/'Dt denotes the material time derivative acting on a tensor, 
say Aij. For present purposes, we can write: 

= — I- h iOkj — U)ik A-kj] (2.5) 

where Uij is the spin: 

Equation (12. 4 p implies that we are neglecting nonlinear elasticity, and are assuming that 
the elastic parts of all displacements are small. Note, however, that we are making no 
such assumption about the plastic displacements. In this Eulerian formulation, neither 
displacements nor strains appear explicitly in the equations of motion, and the velocity field 
Vi{x,t) can describe arbitrarily large and complex motions of material points. 

Consistent with our assumption of constant density po, we assume that the plastic part 
of the rate-of-deformation tensor Dfj, like Sij, is a traceless symmetric tensor, so that plastic 
deformation is volume conserving. This assumption is not necessary for any theoretical 
purpose; it simply reduces the mathematical complexity of the analysis and is an accurate 
approximation in most physical circumstances. When dilation becomes important, it is easy 
to add a volume nonconserving term to Dfj. 



B. STZ Equations of Motion 



The role of the STZ theory is to provide a constitutive relation between the plastic rate 
of deformation tensor Dfj and the deviatoric stress tensor Sij. For illustrative purposes, we 
focus on a common class of situations in which plastic deformation is slow on molecular time 
scales, and in which Dfj has the form 

ToDfj = e-^^/^fijis,e). (2.7) 

Here, s denotes the stress tensor; 9 = ksT and x = ^sTeff are, respectively, the ordinary 
and effective temperatures in energy units; ez is an STZ formation energy; and tq is the 
molecular time scale. A low-temperature expression for fij{s, 6) is shown below in Eq.f l3.30|) . 
In principle, -Dfj also depends on two internal state variables: a dimensionless STZ density 
A, and a tensor that carries orientational memory. Both of these variables must satisfy 
their own equations of motion, and both play central roles in the derivation of the STZ 
equations of motion in Sec JIIIDl In going from such equations to the special form shown 
in Eq. fl2.7p . we have assumed that both A and rriij have reached their equilibrium values 
in times short compared to the time scale for plastic deformation. As will be seen, it is 
the dynamics of rriij that, at low temperatures and for quasi-static deformations, causes 
fij{s, 6) to vanish when the magnitude of the stress |s| is smaller than a dynamic yield stress 
Sy. At somewhat higher temperatures, the m dynamics produces a smooth transition from 
thermally assisted creep to steadily driven flow at about the same yield stress. By dropping 
rriij on the right-hand side of Eq.f l2.7p . we exclude consideration of Bauschinger effects or 
anelastic stress-strain responses below the yield stress, which are contained naturally in more 
general versions of this theory. 

Equation f l2.7p reflects the central premise of the STZ theory that, during plastic de- 
formation, localized, irreversible molecular rearrangements take place sporadically at well 
separated sites. The rearrangements that have been observed directly in numerical simula- 



tions 



fly, 



40| are STZ transitions from one of their orientational states to another. The 
STZ's themselves are ephemeral, noise activated, configurational fluctuations that happen 
to be susceptible to stress-driven, shear transformations. An STZ, when formed, rapidly 
undergoes a shear transition if it is aligned favorably with respect to the stress. Once this 
happens, it cannot transform further in the original direction; but it can transform back- 



wards if the stress is reversed before significant further deformation occurs 



In either 



case, the STZ ultimately disappears into the background of noisy fluctuations. This picture 
of STZ transitions as infrequent events, and the correspondingly long time scales associated 
with plastic deformation, is expressed in Eq.f l2.7p by small values of the Boltzmann-like 
factor exp { — ez/x)- 

Equation fl2.7p must be supplemented by an equation of motion for x- As will be seen in 
Sec JIII Cl this equation is fundamentally a statement about entropy flow in a driven system. 
It has the form 

ez 

where ki and H2 are dimensionless constants. The flrst factor on the right-hand side of 
Eq. (l2.8p . i.e. the product exp (— e^/x) r(s, 6*), is proportional to the rate at which entropy 
is produced by the driving forces. We show in SecJIIII that, at low temperatures, this 
product is proportional to the work rate D^jSij, which contains the factor exp(— e^/x) 
according to Eq.( l2.7p . Thus, the plastic deformation and the effective temperature are 
similarly slow quantities. More generally, r(s, 6) is a non-negative noise strength that must 
be temperature dependent, because the work rate by itself can be negative if, for example, 
thermal fluctuations drive plastic flow in a direction opposite to the stress 

The second factor in the term proportional to ki tells us that x approaches a steady-state 
value Xssil)- Here, q = tq |Dp'| = tq ^JJl/2yDT^ is the magnitude of the plastic strain 
rate in units of the molecular frequency Tq"^. Under normal circumstances, g ^ 1, and 
Xss{0) = Xo is a measure of the disorder induced by slow straining or stirring. On the other 
hand, at very large strain rates, g ~ 1, and Xssil) becomes large. 

The term proportional to /«2 on the right-hand side of Eq.f l2.8p is the rate at which x 
relaxes to 6 in the absence of external driving. This rate contains the factor exp (— e^i/x), 
which determines the frequency of conflgurational fluctuations that couple to ordinary ther- 
mal fluctuations, in rough analogy to the way in which STZ's couple to the external stress. 
In general, we expect the formation energy for such fluctuations to be different from ez- 
The factor p{6) is the attempt frequency for the thermal coupling events. It is a super- 
Arrhenius function of 6 that we presume to vanish below the glass transition, implying that 
the effective temperature is frozen into the system at sufficiently low 9, where aging ceases. 
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FIG. 1: Theoretical stress-strain curves for Vitreloy 1, at 643 K, for different strain rates as shown. 
These curves are in good agreement with data reported by Lu et al. [4l[. The one exception is that, 



for the top-most curve, there is no data at strains beyond the stress peak, presumably because the 
sample failed at that point. 




FIG. 2: Tensile stress for Vitreloy 1 as a function of the scaled strain rate 2?77v7, where t/tv is the 
N^ewtonian viscosity. The data points, with temperatures as indicated, are taken from Lu et al. 



41[. The three solid red curves, from bottom to top, are theoretical predictions for temperatures 



T = 573 K, and 683 i^T respectively. 
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C. Stress-Strain Relations for a Bulk Metallic Glass 



As a first illustration of the STZ theory in operation, consider the deformation 



measurements carried out by Lu et al. 



41| using the bulk metallic glass Vitreloy 1 



(Zr4i,2 Tii3 8 Cui2.5 Niio Be22.5)- In these experiments, a uniform bar of this material was 
subjected to a uniaxial compressive stress, which was measured as a function of strain over 
a wide range of constant strain rates, and over a range of temperatures above the glass 
transition. The STZ analysis of this data is described in . Here, we summarize only 
general features of the comparison between theory and experiment. 

Theoretical stress-strain curves for four different homogeneous strain rates are shown 
in FigHJ This set of curves, and a similar set for different temperatures, are all in good 
quantitative agreement with the experimental data. As seen here, the stress first rises 
elastically, proportional to the strain, while Dfj remains small on the right-hand side of 
Eq.f i2.4p . As X density of STZ's increases according to Eq. fl2.8p . the plastic flow 

becomes dominant, and the stress relaxes to its steady-state value. Almost all of the STZ 
parameters used in plotting these curves were determined from steady-state data, including 
those appearing in the transition-rate formula, shown below in Eq. (13.29 p . and the values 
of the thermal coupling factor p{6) in Eq.f l2.8p . which were obtained from the measured 
Newtonian viscosity. The transient behavior in Fig{Tl i.e. the crossover from elastic to 
plastic response, is determined primarily by the parameter ki in Eq.f l2.8p or, equivalently, 
the dimensionless effective specific heat c'^^-^ defined below, in the text following Eq.f l3.9p . 
The most important feature of these results is that, with a single fixed value of c'^-^-^ of the 
order of unity, the theory naturally reproduces the positions of the stress peaks and the 
rates at which these transients relax toward steady-state. In other words, the STZ theory 
accurately predicts the nonequilibrium dynamics of these systems, including the competition 
between elastic and inelastic mechanisms, over a broad range of experimental conditions. 

In FiglJl we show a comparison between theory and experiment for steady-state stresses at 
different temperatures, as functions of the strain rate multiplied by the Newtonian viscosity. 
These steady-state stresses are the same as those seen in FigH] in the limit of large strain. 
When the strain rate is scaled in this way, all of the data in the viscous limit of small 
stresses and strain rates automatically falls on a single curve with constant slope. The 
important feature of this figure is that the curves cross over from linear viscosity to what is 
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called "super-plasticity" at increasing strain rate and/or increasing temperature. The full 
curves almost, but not quite, collapse onto each other; the crossover occurs at somewhat 
higher stress for higher temperatures in both the theoretical curves and the data. All of the 
temperatures shown in FigJ2] are above the glass temperature; thus the linear viscosity at 
small stresses can be understood as thermally assisted plastic flow. The nonlinear response 
at larger stresses occurs at approximately the low-temperature yield stress, indicating that 
the yielding mechanism described below, following Eq.f l3.27|) . becomes operative in this 
regime. Thus, the quantitative agreement between theory and experiment in Fig|2] is a 
fairly stringent test of a central feature of the STZ theory. 



D. Shear Banding 

One of the most important applications of STZ theory has been in explaining the mech- 
anism of strain localization that leads to the formation of shear bands and shear fracture. 
This failure mechanism is the primary reason why metallic glasses exhibit limited ductility 
and are not widely used as structural materials. In steels undergoing high rates of deforma- 
tion, shear localization apparently is caused by a feedback of some kind between a softening 
mechanism and the heat released during deformation. In our opinion, this mechanism has 
yet to be understood. In metallic glasses, however, localization is common even at relatively 
ow loading rates; and it has long been suspected, and recently demonstrated convincingly 



42| . that the instability leading to localization must be quite different from that in poly- 
crystalline materials. The thermal conductivity of metallic glasses is too high for adiabatic 
heating to account for an instability on the small length scales observed experimentally. 

In the STZ theory as summarized in Eq.f l2.7p . softening occurs due to variations in the 
effective temperature, which must diffuse only very slowly at rates proportional to the local 
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-|47| provides an ideal 



shear rate. In fact, the existence of shear bands in simulations 
virtual laboratory for testing some of the assumptions of the STZ theory. Since the shear rate 
varies by orders of magnitude from inside to outside the shear band, the effective temperature 
also must vary significantly. Under simple shear-loading conditions, this variation provides 
an opportunity to measure the deformability as a function of structure under laboratory- 
scale applied stresses. Shi and Falk j48| have simulated a two-dimensional, low-temperature, 
binary Lennard- Jones system in order to test the Boltzmann-like relation between effective 
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FIG. 3: (a) Simulated strain rates, averaged over increments of 100% total strain, as functions of 
position at various strains. The dashed line is the imposed average strain rate, (b) STZ predictions 
corresponding to the simulation data in (a), (c) Simulated potential energy per atom as a function 
of position at the same total strains shown in (a), (d) STZ predictions for the effective temperature 
as a function of position. The dashed lines in (c) and (d) show the initial values for the potential 
energy and effective temperature. 



temperature and shear rate in Eq. fl2.7p . To do this, they used the "quasi-thermodynamic" 
a^^umption fl that the effec^ve temperature . proportional to the average potential energy 
per atom. Manning, et al. [49| have solved the STZ equations shown in Sec JII Bl above, and 
have found good agreement with the molecular dynamics results as well as some unexpected 
interpretations of them. 



The comparisons between simulation 48| and theory 49| are shown in Figj3l The system 
is a two-dimensional strip subject to simple shear tractions imposed along the upper and 
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lower edges. The theory used a simphfied athermal STZ transition rate R{s), defined below 
in Eg. (13. lip , that rises linearly at stresses appreciably larger than Sy. Figures [3] (a) and (b) 
show, respectively, the simulated and theoretical shear rates, averaged over the length of the 
strip, as functions of position along the transverse direction denoted by Y. As indicated, the 
different curves are snapshots at different total strains ranging up to 800%. FigslH] (c) and 
(d) show the potential energy and effective temperature as functions of position at roughly 
the same sequence of total strains. In accord with the quasi-thermodynamic assumption, 
these sets of functions track each other accurately. 

The quantitative agreement between the simulations and theoretical results shown in 
FigJ31 along with a stability analysis in 49| , reveals that shear banding in these materials is 
a nonlinear, transient instability. The system is initially in a state of uniform shear indicated 
by the dashed horizontal lines at the bottoms of Figsj3] (a) and (b). The dashed curve at 
the bottom of Figj3] (c) is the initial potential energy, whose irregularity was determined 
by the rate at which the sample was quenched from a high temperature. The irregular, 
initial effective temperature in Figj3] (d) was chosen to have approximately the same spatial 
noise spectrum as the simulated potential energy. The band appears only when this spatial 
irregularity has a large enough amplitude, and then only when the strain rate and the initial 
average effective temperature satisfy conditions discussed in [49|. Its position depends on the 
initial noise distribution; but the numerical and theoretical bands behave almost identically. 
They rise rapidly and, for a while, take up almost the entire strain rate, which drops to a 
very small value outside the bands. Both the potential energy and the effective temperature 
saturate inside the bands at limiting values corresponding to Xss(O) = Xo i^i Eq.f l2.8p . At very 
late stages, when the total strain has reached multiples of 100%, the band slowly spreads 
out and collapses, because the small strain rates in the outer regions of the system slowly 
drive x to its steady-state value Xo; and the entire system flows plastically. 

In general, the steady-state effective temperature Xssil) rises rapidly when the dimen- 
sionless strain rate q approaches unity. According to the analysis in 50| , this property of Xss 
causes strongly driven shear bands to collapse, producing very narrow, fracture-like, failure 



zones. Daub and coworkers 



51 



52| have used the STZ theory to describe the dynamics of 



the granular material in an earthquake fault, and have shown that this fracture mechanism 
can account for the sudden stress drops sometimes observed in large seismic events. 
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E. Free-Boundary Problems 



Perhaps the most ambitious goal of the STZ theory is to use it in the full, elasto-plastic 
equations of motion shown in Sec lII At and to predict time-dependent deformations of finite 
systems subject to external tractions. A first, numerically unsophisticated step in this 
direction was made in |53|. More recently, Bouchbinder and coworkers have used STZ 



plasticity in a studies of cavitation instabilities. |54j-l56| The computational problem is 
challenging, partly because including both rapid elastic and slow plastic responses in a single 
numerical procedure is difficult, and especially because this is necessarily a free-boundary 
problem in which the geometry is changing as a function of time. 

Figure m shows recent results by Rycroft and Gibou 57|, in which a necking instability is 
followed all the way to fracture. The system is a two-dimensional strip subject to tractions 
exerted by infiexible, vertically sliding grips on the left- and right-hand sides, which move 
away from each other at a fixed speed. The strip initially has a smooth notch near the center 
of its upper edge. The red and blue regions indicate higher and lower effective temperatures 
respectively. In the top picture, a pair of effectively hotter, i.e. internally disordered, 
shear bands emerges from the notch along the directions of maximum shear stress. These 
bands, like the one shown in Figj3l then broaden into slipping regions. Ultimately, the strip 
separates into two parts, each with a pattern of residual disorder in the places where the 
local plastic fiow was largest. Other results, not shown here, indicate that elastic energy 
initially is stored uniformly throughout the strip. Then, as the necking instability grows, 
this energy fiows to the neck and is dissipated there. 
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FIG. 4: Four snapshots of a necking instability computed by Rycroft and Gibou js^] by solving the 
elasto-plastic equations of motion with STZ plasticity. The red regions are effectively hotter; i.e., 
they have higher effective disorder temperatures, and therefore have undergone more irreversible 
plastic deformation than the bluer regions. 

A detailed description of the numerical procedure used to generate these pictures can 
be found in [57[. Very briefly - The simulation was based on an athermal STZ theory 
of the kind described here in Sec lIIIl with a very simple rate factor R{s) similar to that 
used in the shear-banding analysis described above. The equations of motion were the two- 
dimensional versions of those shown in Sec lII At except that a small viscous term proportional 
to V^f j was added to the right-hand side of Eg. (12.21) in order to damp out elastic oscillations. 
The boundary was tracked using a level-set method. The upper and lower edges were free 
surfaces; there was no surface tension. 

The main limitation of this numerical scheme is that, so far, it has been useful only for 
describing ductile behavior of the kind seen in these pictures. Had the model been brittle, 
or perhaps had it been numerically possible to explore substantially larger pulling speeds 
with the same model parameters, one or more cracks might have started at the notch and 
propagated downward through the system. If this technical limitation can be overcome, we 
should have a powerful tool for studying dynamic fracture. 
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III. THERMODYNAMIC DERIVATION OF THE STZ EQUATIONS 



A. Athermal STZ Model 

For simplicity, in this part of the review, we focus primarily on what we call the "ather- 
mal" limit of plasticity theory. By "athermal," we do not mean strictly zero temperature. 
On the contrary, as stated in SecUl we assume that there is always some thermal or mechan- 
ical noise that sets the time scale tq for rapid, small-scale motions. However, we assume that 
this noise is not strong enough to cause large-scale molecular rearrangements in the absence 
of external forcing. In particular, p{9) = in Eq. fl2.8p . As a result, models of this kind 
have well defined yield stresses, but not linear viscosities, and they do not exhibit thermally 
induced strain recovery. They do describe, for example, irreversible deformation of glasses 
below or near the glass temperature, or the flow of densely packed granular materials sub- 
ject to stresses large enough that they become unjammed. At the end of this Section, we 
discuss briefly how the athermal theory has been supplemented to produce the more general 
equations of motion shown in SecJITl 

Our thermodynamic analysis is based on the assumption that the configurational degrees 
of freedom within a solidlike material are driven out of equilibrium with the heat bath when 
the system is persistently deformed by external forcing, and that they are naturally described 
by an effective temperature under those circumstances. {27] We use the term "configurational 
degrees of freedom" to denote the mechanically stable molecular positions that change slowly 
during irreversible deformation, as opposed to the much faster molecular vibrations about 
the stable configurations. Mathematically, the configurational degrees of freedom specify the 
"inherent structures" 
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60|. Formation of STZ's, or STZ transitions between their internal 
orientational states, are events in which the system moves from one inherent structure to 
another. 



B. First and Second Laws of Thermodynamics for a Plastic Solid 

The preceding discussion implies that an amorphous, solidlike material consists of two, 
weakly coupled subsystems: the slow configurational degrees of freedom on the one hand, 
and the fast kinetic- vibrational degrees of freedom on the other. The fast degrees of freedom 
are strongly coupled to a heat bath so that, together, they and the heat bath constitute a 
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thermal reservoir at temperature 9 = ksT. 

It is useful to start with a microcanonical formulation in which the energy Uc of the 
configurational subsystem is a function of its entropy Sc, its volume V , an elastic shear 
strain e, and a set of internal variables {A} that in Sec Jill Dl will become the number density 
of STZ's and a measure of their average orientation. Throughout this discussion, we choose 
entropies to be dimensionless quantities - logarithms of numbers of states - and express the 
temperatures 6 and x iii units of energy. For simplicity, consider only pure shear deformation 
in, say, the x, y plane, so that the deviatoric stress tensor has components s^x = ~ Syy = s, 
the elastic strain tensor is e^x = ^ ^yy = ^i and the rate of plastic deformation tensor 
is D^^ = — D^y = D^K Let the thermal reservoir have energy Ur and entropy Sr. This 
reservoir has none of its own internal degrees of freedom, and does not support a shear stress. 
According to the definition of temperature, the effective temperature of the configurational 
subsystem is 

which is not necessarily the same as ^ = 81/^/83 r. 
The total energy of this system is 

U"^' = UciSc,e,{^}) + URiSn). (3.2) 
The first law of thermodynamics, 

2VsL>*°* = f7*°*, (3.3) 

says simply that energy is conserved when work is done on the system at the rate 2V s 
Assume, as in Eg. (12.41) . that the total rate of deformation is the sum of elastic and 
plastic parts, i.e. = e + D^''. li V s = {dUc/de)sc,{A}, i-e. if the stress is wholly elastic 
in origin, then the elastic terms cancel out on either side of Eq. (l3.3p . We therefore omit e, as 
well as the constant volume V, as explicit arguments of Uc and Sc- The first law becomes 



R- (3.4) 

Sc 



The fundamental statistical statement of the second law of thermodynamics is that the 
total entropy of an isolated system is a non- decreasing function of time: 

Stot = Sc + Sr>0. (3.5) 
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As argued in 26|, this statement is thermodynamically self-consistent only if the set {A} 
consists of a small number of state variables, each of which is an extensive quantity (or the 
volume average of such a quantity) . Using Eq. (13.41) to evaluate x Sc in Eq. (I3.5p , we find 



W{s,{A}) + {x-0)Sn>O; W{s, {A}) ^ 2V s D^^ -^^ (^] A,. (3. 



6) 



This inequality must be satisfied for arbitrary, independent variations of the A^ and Sr] 
thus each of its component terms must separately be non-negative. We immediately enforce 
{x ~ 0) Sr > by writing 



R 



9Sr = A{x,0) {x-0) = -Q, 



(3.7) 



where A{x, 0) is a non-negative thermal conductivity, and Q is the rate at which heat is 
flowing from the thermal reservoir into the configurational degrees of freedom. 

The inequality VV(s, {A}) > is a form of a Clausius-Duhem inequality that requires a 
non-negative rate of heat production; that is, the rate at which work is done must exceed 
the rate at which energy is stored internally. We will use this inequality in Sec JIII D] to 
deduce features of the STZ equations of motion. Note that we have derived this inequ ality 
from fundamental principles, using an unambiguous statistical definition of the entropy 26| , 
rather than postulating it as an axiomatic form of the second law. The latter strategy is the 
one that is common in the literature. See, for example, the monographs by Lubliner 61| . 
Maugin |62[|, and Nemat-Nasser 



and Gurtin 



64 



63| : or the classic series of studies by by Coleman, Noll, 



65]. 



C. Equation of Motion for the Effective Temperature 

Our first-law equation, Eq. (l3.4p . now has the form 

xSc = W{s,{A}) + Q. (3.8) 

We can use Eq. (13.81) to derive an equation of motion for x by making several observations. 
First, although the STZ's account for all of the coupling beween the applied stress and the 
plastic deformation, they are very rare fluctuations and constitute only a negligibly small 
fraction of the total energy or entropy of the configurational subsystem. Thus Eq. (13.81) is a 
simple statement of energy conservation that can be reduced to 

Vc^^^X-'^VsDP^ + Q. (3.9) 
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where V d^^^ = x {dSc / dx) is the effective heat capacity. 

Second, note that the only relevant rate factor in this athermal system is the work rate 
2V s D^^ itself. So long as there are no thermal fluctuations capable of inducing reverse 
plastic flow, D'P^ must have the same sign as s, and this rate is non-negative. Moreover, in 
the absence of such fluctuations, Q must be proportional to 2V sD^''] the configurational 
system does not move at all without external forcing. 

Third, we know that x must reach some steady-state value during steady shear flow. As 
in Eq. fl2.8p . define the dimensionless strain rate q = tq |-D^'|, and denote the steady-state 
effective temperature by Xssil)- (See [6g] for a detailed discussion of the q dependence of 
Xss-) In the limit g ^ 1 for an athermal amorphous system, Xo = Xss(O) is roughly (perhaps 
exactly) equal to the glass transition temperature, i.e. Xo ~ Tg. In other words, athermal 
systems reach fluctuating steady states of disorder under slow shear. The slower the shear, 
the longer the system takes in real time to reach steady state; but the ultimate value of 
Xss must be independent of q simply for dimensional reasons - there are supposedly no 
competing time scales when g — ?► 0. By definition, the right-hand side of Eq. fl3.9p vanishes 
when X = Xssil)- Therefore, for x not too far from Xssig), we approximate Eq.f l3.9p by 

X 

Xss{q)_ 

Here we see explicitly that the characteristic time scale for x is the same as the time scale 
for plastic deformation, and both are slow because D^'' is proportional to the small density 
of STZ's. 



1 - 



(3.10) 



D. STZ equations of motion 

We turn now to constructing an athermal STZ model based on effective-temperature 
thermodynamics. 

It is easiest and physically most transparent to assume that the STZ's are oriented only 
in the ± directions relative to the stress. In fact, we lose no generality by doing this; the 
tensorial generalizations of the equations are obvious at the end of the analysis. Let the 
number of ± STZ's be N±, and let the total number of molecular sites be N. The master 
equation for the N±. is 

To iV± = R{±s) - R{ts) iV± + r{s) - iV±) , (3.11) 

19 



where R{±s) / tq is the rate factor for STZ transitions between their orientations, and r(s) / tq 
is the corresponding factor for noise driven creation and annihilation of STZ's. The equi- 
hbrium number N^g and the rate factor T{s) will be determined shortly by thermodynamic 
arguments. The internal state variables introduced in Sec JIII Bl are 

, +N_ N.- 



A 



m 



(3.12) 



Here, A is the fractional density of STZ's; and m is their orientational bias which, as men- 



tioned following Eq. (12.7p . becomes the traceless, symmetric tensor rriij in more general ver- 
sions of the theory. According to Eq. 03.111) . the equations of motion for A and m are 



roA = r(Ae,-A); A 



eq 



N 



and 



where 



Tom = 2C{s) r{s 



m 



P t-qA 
1 m ■ — m, 



A 



C{s) 



R{s) + R{-s) ; r(s 



R{s)-R{-s) 
R{s) + R{-s) 



The rate of plastic deformation is 



R{s) N_ - R{-s) = eoAC(s) r(s) -m 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



where vq is a molecular-scale volume that sets the size of the plastic strain increment induced 
by an STZ transition. We expect eg = Nvq/V to be a number of the order of unity. 

Our model of rare, noninteracting STZ's implies that we can write the entropy in the 
form 

Sc{Uc,A,m) =NSo{A) + N AtPim) + Si{Ui), (3.17) 

where Si and Ui are, respectively, the entropy and energy of all the non-STZ degrees of 
freedom; iplm) is the internal entropy associated with STZ's of average orientation m; and, 
for A < 1, 

^o(A) = -A In A + A. (3.18) 
With this assumption, the configurational energy becomes 



Uc{Sc,A,m) = NAez + Ui{Si) 

= N Aez + Ui(Sc - NSo{A) - N AtPim 



(3.19) 
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We now evaluate the partial derivatives of Uc in Eq. (l3.6p . obtaining 

1"0 



N 



W(s, A,m) 



dFz : ^ . di, 

— — To A — i X A m 7- — 
aA dm 



+ 2AC(s) T{s 



m 



where 



Fz{A, m) = ezA-X SoiA) - X A 



dip 1 

svq + Xtt- > 0, 
om 



ipim) — m — — 
am 



(3.20) 



(3.21) 



As before, the three terms in this inequality must separately be non-negative; but the 
argument, especially for the third term, is nontrivial. The term proportional to A is non- 
negative if 



To A oc — 



dFz 
dA 



(3.22) 



or, more generally, if A has a dynamical fixed point at a minimum of the free-energy-like 
function Fz- This minimum occurs at 



A = A 



eq 



ij[m)e 



-ez/x. 



uim] 



exp 



If \ 

ip(m) — m — — 
dm 



(3.23) 



which is consistent with the definition of Agg in Eq.f l3.13]) . 

The internal entropy ipim) is necessarily a positive, symmetric function with a maximum 
at m = 0; therefore, the second term in Eq. (13.201) is automatically non-negative given a 
properly chosen ip{m). 

The last term in Eq. (l3.20p is the most interesting because, unlike the A term, this in- 
equality does not lead to a free-energy minimization law. Nor does it imply normal flow in 



a free-energy landscape as advocated in 



61| or 63| . It is the only one of the three terms in 



Eq. (l3.20ll that depends explicitly on the stress s, which can, in principle, be assigned any 
value independent of A or m. This term can be made to be non-negative for all values of s, 
and for — 1 < m < 1, by choosing 



dip 
dm 



X 



i{m) 



(3.24) 



where ^{m) is the functional inverse of T(s); that is, T{^{m)^ = m. This choice means that 
both s-dependent factors in this product are monotonic functions that vanish at the same 
m-dependent value of s. We can use this second-law constraint in either of two ways. In 



28|, it was assumed that the STZ's were strictly two-state systems with no internal degrees 
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of freedom, and therefore had an Ising-hke entropy. In that case, the rate factor R{s) had 
to be proportional to exp{vos/x)- A more reahstic interpretation is that the STZ's are 
complex, many-body systems with many internal degrees of freedom. The better strategy, 
then, is to choose a physically motivated form of R{s) and to let that determine ■?/'(m) via 
the choice of Eq.^;2M- 

The latter strategy works especially well in the athermal limit that we are considering 
here. The physically most important feature of that limit is that the rearrangement transi- 
tions always go in the direction of the stress; the noise is not strong enough to drive them in 
the opposite direction. This means that ^ and T(s) ~ sign(s). Equation 

then implies that, for -1 < m < 1, dip /dm ^ and, in Eq. dX^ . 

z/(m) ^ z/(0) = exp [tp{0)] (3.25) 

is the number of molecules in an STZ. Interestingly, the athermal choice of R{s), via the 
second law of thermodynamics, recognizes that ergodicity is broken on time scales relevant 
to STZ transitions. Equation fl3.25p implies that any given molecular site has z/(0) different 
ways of being part of an STZ of size z^(0), independent of the average STZ orientation m. 
This can be true only if the STZ is not switching back and forth between its orientations 
during the time over which we are averaging to compute m. 

Now return to Eq. fl3.20p . With our athermal assumption, only the last term in the ex- 
pression for W remains nonzero. The term proportional to A vanishes for slow deformations; 
and the second term vanishes because dip /dm ^ 0. Up to a factor with the dimensions of 
energy, the quantity W is the non-negative rate at which configurational entropy is being 
generated. It was first argued by Pechenik 3, ^ that the noise strength F should be propor- 
tional to this rate of entropy generation per STZ, with the proportionality factor, say vq sq, 
necessarily having the dimensions of energy. Therefore, 



^0 



N 



TAvoSo = ^W^2AC{s) ris)-mvos; T ^ 2C{s) r{s) - m —. (3.26) 



s 



So 



The resulting relation between the STZ production rate, T{s) Neq/ro in Eq. fl3.1ip . and the 
work rate 2 s D^^ was guessed in 2], and has been confirmed by Heggen et al. jc^] in 
the context of conventional flow-defect theories. This identification of the rate at which 
configurational disorder is created with the strength of mechanically generated noise has 
proved to be a very useful concept, as will be seen below in evaluating the yield stress in 
Eq.dS^HD. 
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The equation of motion for m, Eq. fl3.14p . with A ^ 0, becomes 

r n / 5 TTl \ 

Tom^2Cis) Tis) -m 1 . (3.27) 

L J V So / 

Note that both Eqs. fl3.13p for A and fl3.27p for m describe relaxation to steady state that is 
fast compared to that of the effective temperature described by Eq.f l3.10p . The factor D^^ 
on the right-hand side of Eq. fl3.10p contains the small factor A; but no such factor appears 
in Eq. fl3.13p or fl3.27p . Thus, we confirm that the STZ variables A and m are dynamically 
slaved to relatively slow changes in s and x, which is the assumption that we used in deriving 
Eqs.dlTD and (El]). 

Equation ( 13.27P is the usual STZ theory result. There is an exchange of stability at the 
stress s = sq. For \s\ < sq, the dynamically stable, steady-state solution of Eq. (13.271) is the 
jammed state with m ~ ±1; and the rate of deformation D^^ is zero. On the other hand, 
for |s| > So, the stable solution is m = so/s, and 

So" 



ToDP^ ^ eoz/(0)e-"^/^C(s) 



signi^s, 

s 



(3.28) 



Thus Sq = Sy is the dynamic yield stress. 

To complete the derivation, we need to choose the rate factor R{s). One possibility that 
has worked well in several applications is a thermally activated rate of the form: 

A(s)" 



R{s) = Ro{s) exp 
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A(s) = Aoe""/^. (3.29) 



where -Ro(s) is a symmetric function of the stress. The exponential form of the barrier height 
A(s) is the simplest possible expression that vanishes for large positive stress, diverges at 
large negative stress, and introduces only a single new parameter p.. For the metallic glass 
calculations in j^, -Ro(s) = Jl + (s/si)^ with si ~ Sy. Equation f l3.29p is consistent with 



the athermal approximation if ^ ^ Ao. 

Two generalizations of the athermal equations derived above are needed in order to 
recover the fully thermal STZ theory shown in Eqs. (12.70 and (12. 8p . So long as we are dealing 
with an isotropic material, the only directional information in the system is contained in 
the deviatoric stress. We can then assume that the plastic rate-of-deformation tensor is 



proportional to Sij/\s\, where |s| = J {1/2) SijSij. Thus, using Eq. (l3.28p . we find 



/,,(s)^eoKO)^C(|s|) 1-^ ■ (3.30) 
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Equation fl3.10p , the athermal version of Eq. (12 .Sp , becomes 



1- ^ 



(3.31) 



Xss{q)_ 

with q = tq \D^^\. Further generahzing these results to fully thermal situations is straight- 
forward but considerably more complicated. The essential step is to recognize that the 
mechanical noise strength F introduced in Eq. (l3.111) must become the sum of incoherent 
mechanical plus thermal noise strengths, i.e. F — i- F + p{6), where p{6) is the same thermal 
term that we introduced in Eq. (l2.8p to account for relaxation in the absence of mechanical 
driving forces. The explicit form of F can be computed using the same second-law argument 
that led to Eq. (l3.26p . Details can be found in 

IV. CONCLUDING REMARKS 

So far as we know, the STZ theory is the only existing mathematical description of solid- 
like amorphous plasticity that starts with realistic molecular models, and uses the principles 
of nonequilibrium thermodynamics to guide the prediction of observed phenomena. To date, 
those phenomena have included: the transition between linear viscosity and superplasticity 
as a function of temperature and strain rate for bulk metallic glasses; the transient as well as 
steady-state parts of the stress-strain curves for real and numerically simulated glass-forming 
materials; the nature of transient shear-banding instabilities in glassy materials; and even 
a quantitative understanding of the granular shearing instability that produces sharp stress 
drops during major earthquakes. 

The formulation of nonequilibrium thermodynamics that emerged during the development 
of the STZ theory recently has been extended to a study of memory effects in 

thermally cycled glass formers, i.e. the Kovacs effect josl. This thermodynamic point of 
view even has provided an accurate account of a remarkably wide range of experimental 



data for dislocation- mediated plasticity in poly crystalline solids [69|. The similarities and 
differences between dislocations and STZ's are interesting in themselves. Dislocations are 
well defined entities, directly observable by electron microscopy and subject to fairly well 
understood, deterministic equations of motion. It is only when large numbers of interacting 
dislocations are driven by external forces into chaotic motion that thermodynamic concepts 
become relevant to them. 
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In contrast, the STZ's have never enjoyed the visibihty of dislocations. The elementary 
rearrangements presumably associated with flow defects have been known for decades; but, 
for systems in the process of deformation, it never has been possible to identify the defects 
themselves before the events occurred. The thermodynamic theory developed here implies 
that, with perhaps a few special exceptions, such prior identifications are impossible for most 
practical purposes. In the present theory, the sequence of STZ creation, shear transition, and 
annihilation is a noise activated process, more nearly akin to nucleation of a critical droplet 
in a supercooled vapor than, for example, to the creation of a dislocation at a Prank-Read 
source. We should no more expect to be able to look at a deforming amorphous material 
and predict where the next STZ event will occur than we should expect to be able to predict 
where the next droplet will form in the vapor. Nor should we worry that the stochastic 
nature of STZ plasticity unnecessarily limits the predictive power of the theory. 

The one important case where a deterministic, dynamical theory of amorphous deforma- 
tion should make sense is in the athermal quasistatic (AQS) limit. At zero temperature, 
using numerical simulation, we might be able to strain an amorphous system gradually, and 
predict where the next rearrangement will occur by looking at nearby saddle points in the 
energy landscape. Once the system has crossed a saddle point, however, we cannot predict 
where the next such event will occur unless we stop straining the system and let it relax into 
its nearest energy minimum before resuming the deformation. This is the AQS numerical 
procedure, which often produces system-spanning, avalanche-like events and size-dependent 
noise spectra. There are many real systems that do behave like this, for example, granular 
materials, foams, or colloidal suspensions sheared so slowly that the mechanical noise gener- 
ated by one event has died out before the next event occurs. These are not what we would 
call "normal" plastic materials; they cannot be described by local constitutive laws like the 
ones discussed here. But it may be interesting to locate the boundary between normal and 
AQS systems, and thus to understand the hmitations of the STZ theory. 

In our opinion, however, it will be more interesting to use the theoretical tools developed 
here to explore "normal" plasticity in broader contexts - in particular, to study a variety 
of dense, complex fluids and biological materials. We need to understand the relations 
between the STZ and SGR theories, and perhaps learn how to combine the strengths of the 
two approaches. We have some new tools for understanding nonequilibrium phenomena; 
and we are optimistic that these tools will lead us to new discoveries. 
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